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ABSTRACT 


Arctic  outbreaks  over  the  Canadian  Western  Plains  during  the 
late  spring  period  frequently  take  the  form  of  a  cold  east-northeasterly 
flow  over  a  warmer,  sloping  surface.  A  mesoscale  numerical  model  is 
developed  in  an  attempt  to  simulate  such  circulations.  Following 
Lavoie  (1972)  the  atmospheric  structure  of  the  cold  air  mass  is  repre¬ 
sented  by  three  layers:  a  constant  flux  layer  in  contact  with  the 
earth's  surface,  a  well  mixed  planetary  boundary  layer  capped  by  an 
inversion,  and  a  deep  stratum  of  overlying  stable  air.  Averaging  the 
set  of  governing  primitive  equations  through  the  depth  of  the  mixed 
layer  yields  predictive  equations  for  the  horizontal  wind  components, 
potential  temperature,  specific  humidity,  and  the  height  of  the 
inversion.  Time  dependent  calculations  are  limited  to  this  layer  by 
parameterizing  the  interactions  between  the  mixed  layer  and  both  the 
underlying  and  overlying  layers.  Precipitation  from  limited  con¬ 
vective  clouds,  and  latent  heat  within  the  layer  are  included  in 
terms  of  mesoscale  variables. 

A  47-6  km  by  A7.6  km  grid  mesh  of  1  3&9  points  covering  the 
Canadian  Prairie  Provinces  is  used  to  represent  the  variables  in  the 
mesoscale.  The  governing  equations  are  solved  numerically  with 
terrain  influences,  surface  roughness,  temperature  variations,  and 
moisture  fluxes  allowed  to  perturb  the  mixed  layer  from  its  initial 
conditions  until  resultant  mesoscale  boundary  layer  weather  patterns 
evolve . 
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The  mean  spring  topographic  precipitation  pattern  was 
successfully  reproduced  by  the  simulated  late  spring  upslope  flow 
with  limited  convective  precipitation.  This  indicates  that  late 
spring  upslope  weather  is  a  contributing  factor  to  the  upward  swing 
in  precipitation  amounts  over  the  Canadian  Western  Plains  during  the 
spring  period.  Mesoscale  planetary  bounda ry- 1 ayer  weather  patterns 
appear  to  exert  a  dominant  control  over  the  location  and  intensity 
of  perturbations  in  the  spring  precipitation  pattern.  The  elimina¬ 
tion  of  surface  heating  significantly  reduced  the  area  and  intensity 
of  the  precipitation.  A  case  study  based  on  observed  initial  con¬ 
ditions  showed  that  the  model  could  reproduce  a  persistent  limited 
convective  precipitation  pattern  maintained  by  upslope  flow  and  that 
a  low-level  trough  exerts  a  marked  influence  on  the  location  and  the 
intensity  of  the  precipitation. 
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CHAPTER  1 


INTRODUCTION 

1  . 1  Preliminary  Comments 

The  numerical  experiment  described  in  this  thesis  represents 
an  initial  attempt  to  simulate  quantitatively  the  temporal  and  spatial 
evolution  of  mesoscale  boundary- layer  weather  patterns  resulting  from 
the  interaction  of  an  east-northeasterly  circulation  of  arctic  air 
with  the  warmer,  sloping  Canadian  Western  Plains  during  the  late  spring 
period. 

For  Canada's  Prairie  Provinces,  spring  may  be  defined  as 
the  period  March  through  May  (Longley,  1972).  However,  spring  is  a 
complicated  season  showing  significant  variations  in  both  time  and 
space.  Snow  generally  covers  the  entire  region  in  early  March;  the 
snow  line  moves  progressively  northeastward  pushing  north  of  the 
Prairie  Provinces  by  late  May  (McKay  and  Thompson,  1968).  Climatic 
records  indicate  that  spring  is  a  time  of  increasing  cloudiness, 
precipitation  and  wind  speeds  (Longley,  1972).  Fig.  1  shows  that  it 
is  also  characterized  by  a  mean  east-northeasterly  surface  flow  over 
all  but  the  western  areas  (Bryson,  1966). 
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Streamlines  of  the  surface  resultant  winds  in  May  (adapted  from  Bryson,  1966). 
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Late  spring,  usually  late  April  or  May,  is  the  period  when 
the  northeast  portion  of  the  Prairie  Provinces  (the  Canadian  Shield 
portion)  experiences  rapid  snow  melt,  and  the  Western  Plains  are  snow 
free.  Although  the  'wet  land  surface  period1  on  the  plains  may  be 
quite  short,  and  little  or  no  vegetative  growth  is  underway  to  supply 
moisture  through  transpiration,  some  water-vapour  sources  are  provided 
by  the  numerous  ponds  and  sloughs  which  have  been  filled  by  the  snow¬ 
melt  runoff.  These  small  water  bodies  are  usually  very  shallow,  and 
they  warm  rapidly  along  with  the  snow-free  land  surface.  A  major 
moisture  source  is  provided  by  snow  melt  on  the  Canadian  Shield. 

Arctic  air  sweeping  down  over  the  Prairie  Provinces  from  the 
east-northeast  encounters  a  relatively  warm  surface.  This  readily 
encourages  evaporation  as  well  as  vertical  fluxes  of  momentum  and 
heat.  Thus,  spring's  progression  across  this  region,  interrupted  by 
arctic  outbreaks,  provides  a  complicated  yet  fascinating  period  for 
study . 

1 . 2  Physical  Geography 

For  the  purpose  of  this  study  the  Western  Plains  is  defined 
as  the  region  between  the  Canadian  Shield  and  the  Rockies.  The  eastern 
boundary  extends  from  the  southeast  corner  of  Manitoba  through  Lake 
Winnipeg,  and  its  northern  limit  is  roughly  along  a  line  running  from 
northern  Lake  Winnipeg  to  Fort  McMurray  then  ear tward  to  the  British 
Col umb i a-Al berta  border  between  the  Hay  and  Peace  Rivers  (Fig.  2). 

The  most  general  feature  of  the  relief  of  the  Prairie 
Provinces  is  the  westward  upslope  from  the  Hudson  Bay  lowlands  in 
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Fig.  2.  Topographic  map  of  the  Prairie  Provinces  (after  Longley,  1972). 
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the  east  to  the  Rockies  in  the  west  (Fig.  2).  Within  the  plains 
region  the  topography  is  by  no  means  uniform.  Many  lakes,  ponds  and 
sloughs  occupy  depressions  in  the  landscape.  Upland  areas  are  found 
throughout  the  region.  The  most  notable  is  the  Manitoba  escarpment 
including  the  Pasquia  (2700  ft),  Porcupine  (2600  ft),  Duck  (2700  ft), 
Riding  (2300  ft),  and  Turtle  (2500  ft)  tablelands.  The  Cypress  Hills 
in  southeastern  Alberta  and  southwestern  Saskatchewan  rise  to  *4800  ft 
or  over  2000  ft  above  the  surrounding  plains.  Alberta  has  a  number 
of  other  upland  areas  such  as  the  Swan  Hills  (3800  ft)  south  of  Lesser 
Slave  Lake  and  the  Birch  Mountains  (2800  ft)  northwest  of  Fort 
McMurray . 

1 . 3  Late  Spring  Upslope  Weather 

Late  spring  on  the  Western  Plains  is  often  dulled  by  extensive 
strata  of  limited  convective  clouds  (stratocumul us) .  Cool,  moderately 
moist  winds  gust  out  of  the  east-northeast  and  the  ground  is  dampened 
by  light  drizzle  showers,  not  infrequently  mixed  with  snow.  This 
type  of  weather  is  not  only  unpleasant  for  those  trying  to  extend  the 
short  summer  recreation  period,  but  also  has  severe  economic  implica¬ 
tions.  In  an  area  where  agriculture  is  still  the  major  industry, 
the  grey  clouds,  blocking  out  the  sun,  retard  surface  drying  and  the 
precipitation  further  moistens  the  earth.  This  tends  to  keep  the 
farmers  off  the  land  during  the  critical  spring  mowing  period.  The 
extensive  low  cloud  also  causes  concern  to  those  whose  interests  are 


in  aviation. 
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Climatological  charts  (Longley,  1972)  indicate  that  spring 
is  characterized  by  increasing  cloudiness  and  precipitation.  Most 
certainly  one  of  the  contributing  factors  is  the  occurrence  of  late 
spring  upslope  weather.  The  term  upslope  generally  refers  to  a 
persistent,  shallow,  east-northeasterly  circulation  associated  with 
outbreaks  of  arctic  air  moving  up  the  gently  rising  slope  of  the 
Western  Plains.  In  late  spring,  this  is  normally  a  relatively  unstable 
situation.  The  cold  air  mass  moving  south-westward  up  the  sloping 
plains  is  heated  from  below  by  contact  with  the  snow-free  ground. 

This  often  produces  turbulent  low-level  winds,  and  extensive  decks  of 
limited  convective  clouds  (st ratocumu 1  us)  accompanied  by  light  drizzle 
showers  or  snowf 1 urr ies .  Limited  convective  condensation  forms 
generally  imply  that  the  atmosphere  is  characterized  by  an  unstable 
stratification  of  limited  depth  which  is  topped  by  a  marked  stable 
stratum  (Willet  and  Sanders,  1959)*  The  air  is  moistened  from  below 
by:  (l)  the  numerous  lakes  making  up  over  half  the  area  of  a  thawing 

Canadian  Shield,  including  lakes  Winnipeg,  Manitoba,  and  Wi nn i pegos i s , 
which  lie  to  the  east  of  the  Western  Plains  and  (2)  the  large  number 
of  ponds  and  sloughs  which  dot  the  plains.  The  unstable  air  flow, 
readily  available  moisture,  and  the  sloping  terrain  must  certainly 
contribute  to  the  increasing  cloudiness  and  precipitation  during  this 
season  of  the  year. 

The  synoptic  situation  which  produces  a  persistent  east- 
northeasterly  flow  across  the  Western  Plains  usually  consists  of  an 
arctic  high  moving  from  the  Great  Bear  Lake  area  towards  Lake  Superior 
in  conjunction  with  a  low  (especially  a  cold  low)  moving  across  the 


7 


northern  United  States. 

The  streamlines  of  the  mean  surface  winds  for  May  (Fig.  l) 
show  an  east-northeasterly  flow  over  all  but  the  western  areas.  The 
maps  for  March  and  April  show  a  similar  flow  pattern  (Bryson,  1 966 ) . 
These  are  mean  maps.  The  day-to-day  situation  often  shows  the  arctic 
front  in  a  trough  just  east  of  the  Rocky  Mountains,  subjecting  the 
entire  Western  Plains  to  an  east-northeasterly  flow. 

Although  the  Western  Plains  are  often  said  to  be  in  the  rain 
shadow  of  the  Rocky  Mountains,  actually  the  Rockies,  and  the  sloping 
terrain  from  the  Hudson  Bay  lowlands  to  the  foothills  of  the  Rockies, 
increase  rather  than  decrease  precipitation  over  the  plains  during 
the  spring  period.  This  has  been  shown  convincingly,  at  least  for 
Alberta  and  southwestern  Saskatchewan  by  Reinelt  (1969).  It  is 
readily  seen  from  the  mean  spring  precipitation  chart  (Longley,  1972) 
that  the  foothills  and  upland  areas  receive  appreciably  more  precipita¬ 
tion  than  the  surrounding  lowland  areas.  This  lends  support  to  the 
proposition  that  unstable  upslope  flow  is  a  contributing  factor  to 
the  increased  spring  precipitation. 

1 . k  Outline  of  the  Study 

This  study  uses  a  three  layer  numerical  model  to  simulate 
the  mesoscale  boundary  layer  weather  patterns  induced  by  a  moderately 
moist,  cold,  east-northeasterly  flow  moving  up  the  relatively  warm 
sloping  Western  Plains.  Time-dependent  calculations  are  limited  to 
a  single  mixed  layer  by  parameterizing  the  interactions  of  this  layer 
with  the  underlying  and  overlying  layers.  The  procedure  is  based  on 
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the  hypothesis  (after  Lavoie,  1972,  197*0  that  a  single  well-mixed 
layer  may  be  used  to  represent  many  of  the  significant  aspects  of 
mesoscale  boundary  layer  weather  disturbances,  and  that  these  meso- 
scale  boundary  layer  disturbances  provide  effective  control  over  the 
location  and  intensity  of  convective  activity. 

The  model's  ability  to  simulate  'typical'  late  spring  upslope 
weather  was  tested  initially.  The  relative  importance  of  surface 
heating  compared  to  friction  and  terrain  relief  in  the  production  of 
limited  convective  precipitation  was  evaluated  by  running  the  model 
while  neglecting  vertical  eddy  heat  fluxes  under  otherwise  identical 
meteorological  conditions.  A  case  study  was  used  to  determine  the 
model's  ability  to  handle  situations  based  on  observed  initial  con- 
d i t i ons . 

The  details  of  this  study  are  presented  in  the  following  four 
chapters.  Chapter  2  outlines  the  governing  equations,  together  with 
details  of  the  model  atmosphere.  Chapter  3  contains  a  brief  description 
of  the  numerical  technique  used  in  this  study,  along  with  a  step-by- 
step  outline  of  the  numerical  procedure.  The  experimental  design, 
together  with  the  various  results,  are  explained  in  Chapter  4.  A 
brief  summary  and  conclusions  make  up  Chapter  5. 


CHAPTER  2 


THE  MODEL 


2. 1  Treatment  of  the  Planetary  Boundary  Layer 

Following  Deardorff  (1972),  the  planetary  boundary  layer 
(PBL)  may  be  defined  as  the  region  of  the  atmosphere  adjacent  to  the 
earth's  surface  where  small-scale  turbulence  is  induced  by  thermal 
convection  and  wind  shear.  Clark  (1970)  and  Deardorff  (1972)  have 
suggested  that  there  are  two  possible  approaches  available  for  the 
treatment  of  the  PBL.  One  is  to  resolve  the  vertical  structure  of 
the  lower  atmosphere  explicitly  by  placing  several  layers  in  the 
lowest  1  to  3  km.  This  usually  involves  (1)  pa rameter i zat i on  of 
fluxes  at  the  surface  and  (2)  modelling  the  turbulent  diffusion 
of  properties  above  Z$,  the  height  of  the  lowest  model  level.  A 
second  approach  is  to  parameterize  all  aspects  of  the  PBL  so  that 
its  structure  does  not  explicitly  come  into  the  computations;  Z$ 
is  set  at  a  much  greater  height,  say  1  km  and  part  (2)  is  by-passed. 

Frequently  the  PBL  is  divided  into  a  constant-flux  (Prandtl) 
layer  in  which  the  fluxes  of  heat,  momentum  and  water  vapour  are 
considered  constant  (i.e.,  having  the  same  magnitude  as  at  the  surface) 
reaching  some  10  to  80  m,  and  a  second  higher  layer  in  which  the 
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fluxes  generally  decrease  with  height.  Under  the  conditions  described 
in  Chapter  1,  a  well-mixed  atmosphere  is  often  present  in  the  lowest 
1  to  2  km  of  the  cold  air  mass.  Mixing  is  induced  both  thermally  as 
a  result  of  surface  heating  along  the  trajectory  of  the  arctic  air, 
and  mechanically  by  moderate  low-level  winds.  Thus,  the  analysis  of 
the  PBL  may  be  simplified  by  assuming  that  overlying  the  constant-flux 
layer  is  a  convective  layer  which  is  well-mixed.  It  is  assumed  that 
the  mixed  layer  is  capped  by  an  inversion.  This  upper  stable  layer 
completes  the  atmospheric  structure. 

2 .2  Thermal  Structure  Assumed  in  the  Model  Atmosphere 

The  atmosphere  is  assumed  to  consist  of  three  layers  defined 
by  distinctive  lapse  rates  of  potential  temperature  (Fig.  3).  The 
vertical  thermal  structure  includes: 

(a)  a  surface  layer  (Zq  to  Z^)  characterized  by  a  super- 
adiabatic  lapse  rate  giving  rise  to  a  condition  of  upward  heat  flux. 

The  upper  limit  is  set  at  10  m,  the  anemometer  level. 

(b)  a  mixed  layer  (Z^  to  h).  in  which  the  horizontal  wind 
velocity  V,  the  potential  temperature  0,  and  the  specific  humidity  q 
are  assumed  homogeneous  in  the  vertical.  Its  upper  limit  is  defined 
by  a  first-order  discontinuity  in  potential  temperature.  The  depth 
of  this  layer  varies  locally  in  space  and  time  as  a  result  of  several 
forcing  factors  including:  (1)  varying  surface  heat  fluxes,  (2)  varying 
surface  momentum  fluxes,  (3)  terrain  slope,  and  (4)  the  release  of 
latent  heat  from  precipitation. 
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Fig.  3. 


Thermal  structure  assumed  in  the  model  atmosphere 
(after  Lavoie,  1972,  197*0  . 
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(c)  an  upper  layer  (Z>h)  consisting  of  a  deep  stable  stratum 
with  a  constant  lapse  rate  of  potential  temperature. 

2 . 3  Mesoscale  Phenomena 

The  model  will  be  applied  to  weather  systems  whose  character¬ 
istic  dimensions  are  of  the  order  of  tens  of  kilometers  in  the  hori¬ 
zontal  and  less  than  ten  kilometers  in  the  vertical.  These  phenomena 
fall  within  the  scale  of  atmospheric  motion  usually  referred  to  as 
mesoscale.  Ogura  (1963)  points  out  that  it  is  customary  to  assume 
hydrostatic  balance  for  motions  whose  horizontal  scale  is  much  larger 
than  the  vertical  scale.  The  assumption  that  the  pressure  field  is 
determined  by  the  density  distribution  alone,  without  reference  to 
the  field  of  motion,  rules  out  the  possibility  of  simulating  individual 
convective  clouds.  However,  it  is  hypothesized  (after  Lavoie,  1972, 
197*0  that  the  location  and  intensity  of  convective  regions  will  be 
determined  by  the  larger  mesoscale  weather  patterns. 

2 . k  Governing  Equation 

The  basic  equations  applied  to  a  fluid  parcel  in  the  mixed 
layer  (Lavoie,  1972,  197*0  are 

Horizontal  momentum  equation 

=  -  k  x  f V  -  aVp  +  a  ~  (l) 

First  and  second  laws  of  thermodynamics 
d0  _  a6  9Q  aL0M " 

dt  "  ’  C  T  9Z  C  T  K  ' 

P  P 


Conservation  of  mass 


J_  da 
a  dt 


V-V 


Conservation  of  water  substance 


= 

dt 


3E  m- 
'  a3Z  "  aM 


Definition  of  hydrostatic  balance 


3p 

01  3Z  =  "  9 


(3) 


W 


(5) 


Definition  of  potential  temperature 


A  a  k  /  k-i 

0  =  R  P*/P 


Equation  of  state 
pa  =  RT 


(6) 


(7) 


where  the  coordinate  axes  are  chosen  with  the  positive  Z  axis 
upward,  the  positive  x  axis  to  the  east-northeast,  and  the  positive 
y  axis  to  the  north-northwest,  and  where 

t  is  the  vertical  flux  of  horizontal  momentum 

Q  is  the  vertical  flux  of  sensible  heat 

is  the  net  condensation  rate  per  unit  volume 
E  is  the  vertical  flux  of  water  vapour 

p_,,  is  a  reference  pressure  of  1000  mb 

p  is  pressure 

V  is  the  horizontal  velocity  vector 


*■ 
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q  is  specific  humidity 

g  is  the  acceleration  due  to  gravity 
a  is  the  specific  volume  of  dry  air 

T  is  temperature 

W  is  vertical  velocity 

0  is  potential  temperature 

t  i  s  t i me 

f  is  the  Coriolis  parameter 

L  is  the  latent  heat  of  evaporation 

R  is  the  gas  constant  for  dry  air 

0^  is  the  specific  heat  of  dry  air  at  constant  pressure 

k  =  R/C 

P 


The  behaviour  of  the  mixed  layer  as  a  whole  may  be  described 
by  integrating  from  to  h,  using  the  assumption  of  vertical  homo¬ 
geneity  of  horizontal  wind  velocity,  potential  temperature  and 
specific  humidity.  Thus,  the  parcel  derivatives  of  these  quantities 
are  independent  of  height. 

(a)  Horizontal  momentum  equation 
Averaging  (1)  through  the  mixed  layer  yields 


dV_ 
d  t 


k  x  fv 


i  r 

aVpdZ  -  a(x. 


(8) 


Assume  that  the  momentun  flux  through  the  inversion  surface  can  be 
neglected.  The  flux  through  the  lower  boundary  can  be  parameteri zed 
by  means  of  the  bulk  aerodynamic  method  (Priestley,  1959)* 


T 


(9) 


, 


where  C-,  is  the  draq  coefficient  and  p  is  the  density  at  level  Z  . 

D  3  s  1  s 

The  horizontal  pressure  gradient  force  can  be  calculated  in 
terms  of  0  by  eliminating  a  between  (5)  and  (6). 


(g/Cp) 


e ' 


do) 


Assuming  that  pressure  is  known  at  some  level  H  within  the  upper 
stable  layer,  integrate  (10)  from  H  down  to  an  arbitrary  level  Z 


within  the  mixed  layer. 

/ 

k  q  k 

PH  ■  ~  P* 


H  -  h  „  6H 
6H  -  6h  "  6h 


Z  -  h 


1/k 


00 


where  0^  is  the  potential  temperature  just  above  the  inversion  base, 
0H  1  and  0Z  =  0. 

Taking  the  horizontal  gradient  of  p^  multiplying  by  -a  and  approx¬ 
imating  Jin  (0^/0^)  by  the  first  two  terms  of  its  Taylor  series 

expansion  (since  0U/0U  is  approximately  one)  yields 

H  h 


-  aVp. 


(«VP)H  +  g 

0H  H  0H 


2g(H  -  h)  ( 

0..  +  0.  ' 


9h  +  9h 


+  ?e 
6 


(12) 


Because  the  final  term  in  (12)  is  the  only  one  that  depends  upon 

height,  the  vertical  integration  of  (8)  is  readily  accomplished. 

Letting  F  =  —  (aVp)u  represent  the  pressure  gradient  at  level  H 
0 1 1  H 


and  0  =  (0^  +  6|_|)/2  represent  the  mean  potential  temperature  of  the 
layer  (H  -  h)  and  using  (12)  to  re-write  (8)  gives 


dV_ 

dt 


-  k  x  f V  -  fu  +  g  ~  h^  V0  +  ■—  (0  -  0)  Vh 
H  0 . 


(h  -  Z  ) 

+  9  20—— *  V0  - 


I  -> 


'D 


(h  -  Zs) 


(13) 


In  order  to  remove  the  critical  dependence  of  the  second 
and  third  terms  of  (13)  upon  the  height  of  the  undisturbed  level  H, 
the  geostrophic  wind  shear  in  the  upper  stable  layer  is  introduced 
in  terms  of  the  thermal  gradient.  The  thermal  wind  equation  is 


g 


0 


V0 


fk  x 


3V 
_ 9 

8Z 


(14) 


where  \l ^  is  the  geostrophic  wind  in  the  upper  layer  and  ip  is  the 
rotated  shear  vector.  If  can  be  considered  essentially  constant 
in  the  upper  layer,  then  applying  the  definition  of  F  at  the  initial 
height  of  the  inversion 

=  -  (H  -  h.)  f*  (15) 

n  1  1 

or 

Fu  +  (H  -  h)  f$  =  F.  +  (h.  -  h)  f^  (16) 

H  1  l 


Substituting  (16)  into  (13)  eliminates  reference  to  level  H,  except 

for  -2-  (0  -  0)Vh  which  represents  the  restoring  force  due  to  the 
^h 

deformation  of  the  inversion  base.  In  order  to  retain  the  influence 
of  the  stable  layer  aloft,  we  take  the  effective  level  H  to  be 
the  height  of  the  inversion  base  h  plus  one  half  the  difference 


* 
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between  the  initial  height  h.  of  the  inversion  and  the  maximum 
altitude  reached  by  the  inversion  surface  as  a  result  of  the  various 
topographical  influences.  That  is,  for  the  purpose  of  defining  0, 
let 


H  ,  h+l(hm-h.) 


(17) 


where  h  is  the  maxi  mum -d i sturbed  height  of  the  inversion  surface, 
m 

Then 

5  ’  (Bh  +  0H)/2  "  6h  +¥  (hm  -  hi>  (l8) 

where  r  is  the  vertical  gradient  of  0  above  the  inversion.  Using  (14), 
(16),  and  (18)  we  may  write  (13)  as 


3V 

at 


-  V*VV  -  k  x  fV  -  F.  -  (h.  -  h) 

i  i 


+ 


/ 


f  (h- 


Zs> 


V0  - 


TiT^ 


D 


(19) 


The  time-dependent  variables  are  V,  h  and  0. 


(b)  Potential  temperature  equation 

Averaging  (2)  through  the  mixed  layer  and  assuming  (— )  = 


yields 


|i  .  v-ve  -  ^  -  Qs) 


0a 


C  T 
P 


+ 


C  (h  -  Z  )T 
P  s 


(20) 
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where  a  and  T  are  the  mean  values  for  the  layer.  is  the  heat  flux 

through  the  inversion  (neglected).  is  the  heat  flux  at  the  lower 

boundary  (Priestley,  1959) 

Qs  ■  VsCH^(9o  ‘  e)  (2I) 

where  Cu  is  the  sensible  heat  transfer  coefficient  and  0  is  the 
potential  temperature  at  ground  level.  M  is  the  total  rate  of  pro¬ 
duction  of  particles  of  precipitation  size  in  terms  of  mass  of  water 
per  unit  area  by  limited  convective  clouds  in  the  mixed  layer  or, 
alternatively,  the  local  precipitation  rate.  The  pertinent  parameters 
are  taken  to  be  (1)  the  depth  of  the  cloud  given  by  the  difference 
between  the  mean  cloud-top  height,  that  is,  the  height  of  the  inversion 
base  and  the  lifting  condensation  level,  or  the  level  to  which  un¬ 
saturated  air  would  have  to  be  raised  dry  ad i abat i ca 1 1 y  to  produce 
condensation  (cloud  base)  and  (2)  the  mesoscale  field  of  vertical 
velocities  near  cloud  base.  The  precipitation  rate  may  be  represented 

by 

M  =  (cW  +  l)N  mm  hr"1  (22) 

Id 

where  N  =  a(D  -  k)2  is  a  representation  of  the  precipitation  from 
limited  convective  clouds  in  the  absence  of  mesoscale  vertical 
velocities  (Fig.  4) .  The  depth  of  the  cloud  D  must  exceed  k  kilometers 
before  precipitation  is  generated  and  the  precipitation  rate  increases 
more  and  more  rapidly  with  increasing  cloud  depth.  Observations  of 
mid-latitude  layer  clouds  by  Mason  and  Howorth  (1952)  and,  more 
recently,  by  Stewart  (1984),  indicate  that  a  cloud  must  be  400  to  500  m 


Cloud  depth  (km) 
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Fig.  b. 


Relationship  between  precipitation  rate 
and  cloud  depth. 
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thick  for  precipitation  to  reach  the  ground.  Therefore,  k  was  set 

at  450  m.  Parameter  a  is  a  measure  of  the  efficiency  of  production 

of  precipitation  within  the  cloud;  it  was  arbitrarily  set  to  1.0. 

WrD  (cm  sec  *)  is  the  mesoscale  vertical  velocity  at  cloud  base.  The 
Id 

parameter  e  was  set  equal  to  0.5  cm-1  sec  so  that  mesoscale  subsidence 
of  2  cm  sec-1  or  more  is  sufficient  to  prevent  precipitation. 

The  prognostic  equation  for  0  becomes 


90 

9t 


-  v-ve  + 


CH 

(H  -  Zs) 


+  LQa 

C  T 
P 


M 

(h  -  Z  ) 
s 


(23) 


(c)  Mass  continuity  equation 

Averaging  (3)  through  the  depth  of  the  mixed  layer  yields 


J_  — 

a  dt 


W  -  W 
->  h  s 


(24) 


The  most  significant  contribution  to  the  change  in  a  of  an  air  parcel 
is  due  to  vertical  motion 


1  ~  £  w 
a  dt  c2 


(25) 


where  c  is  the  isothermal  speed  of  sound.  Averaging  (25)  through  the 
mixed  layer  gives 

“  ^2  <wh  +  ws> 


a  dt 


(26) 


Using  (25)  and  (26) 


(l  +  g(h  -  Zs)/2c2) 


(h  -  Z, 


X 


(h  -  Z  ) 

9  - 

2c2 


(27) 


where  is  the  vertical  velocity  through  which  can  be  taken  to  be 
entirely  due  to  the  sloping  terrain 

W  =  V-VZ  (28) 

so 


Assuming  that  the  fluid  surface  at  the  top  of  the  mixed  layer  can  be 
considered  to  be  a  material  surface,  then 

37  =  wh  <29> 


In  order  to  prevent  the  development  of  a  superad iabat i c  layer  at  the 
top  of  the  mixed  layer  when  the  inversion  is  wiped  out  by  heating 
from  below  equation  (29)  is  modified  as  follows 


dh 

dt 


wh  + 


0=0, 


where  r  is  the  lapse  rate  of  potential  temperature  in  the  upper  layer 

1  90 

and  -  97  's  the  rate  of  increase  in  the  height  of  the  top  of  the 
mixed  layer  that  will  maintain  a  first  order  discontinuity  in  0  at 
level  h  (included  only  when  0=0^). 


The  prognostic  equation  for  h  becomes 


3h 

dt 


-  V*  Vh  +  W.  + 
h 


I  $1 

r  9t 


(30) 


(d)  Moisture  continuity  equation 
Vertical  integration  of  {k)  yields 


dq  _  a  /r  _  c  \  aM 

dt  "  (h  -  Z  )  n  ’  s'  "  (h  -  Z  ) 

s  s 


(30 


Parameterizing  the  moisture  flux  through  the  lower  boundary 


E 

s 


pscEl^l (q0  -  q) 


(32) 


and  retaining  the  moisture  flux  through  the  top  of  the  mixed  layer 
because  inversions  generally  display  dry  air  aloft 


Eh  =  ph3lwhl (q  "  qh} 


(33) 


where  $  is  an  arbitrary  parameter  (0.5)  and  q^  is  the  specific 
humidity  just  above  the  inversion  base. 

The  prognostic  equation  for  q  becomes 


9q_  = 

St 


-  V-Vq 


cE|v| 

+  (h  -  Zs) 


Blwh'  ,  . 

'  (h  -  Z  )  (q  '  V 

S 


aM 

tf-tt 


where  M  is  given  by  (22)  and  by  (27). 
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(e)  The  complete  set  of  prognostic  equations 


3V 

at 


-V 


V-VV  -  k 


(h.  -  h)  f\p 

i 


z  )ve 

s 


->  I  -> 


(h^zp 


(19) 


30 

at 


-  v*ve  + 


'H 


(h 


-v 


M 

(h  -  Zj 


(23) 


(30) 


ia  = 

at 


-  V • Vq  + 


cEl*| 


BKI  ,  . 

'  (h  -  z  )  (q  '  V 
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aM 

(Fwp 


(34) 
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In  employing  the  above  set  of  prognostic  equations,  the 
hypotheses  and  assumptions  upon  which  they  are  based  must  be  kept  in 
mind.  The  guiding  hypotheses  for  the  simulation  are:  (a)  that  the 
model  atmosphere  is  a  reasonable  representation  of  the  actual  atmos¬ 
pheric  structure,  (b)  that  a  single  well-mixed  layer  may  be  used  to 
reproduce  planetary  boundary-layer  weather  patterns,  and  (c)  that 
mesoscale  PBL  weather  patterns  control  the  location  and  intensity  of 
convective  activity.  The  acceptance  of  these  hypotheses  allows  some 
simplifications  of  nature.  The  simplifications  include:  (a)  the 
assumption  of  vertical  homogeneity  of  V,  0  and  q,  (b)  the  use  of  the 
hydrostatic  approximation,  (c)  the  mesoscale  parameterization  of 
convective  activity  and  (d)  the  limitation  of  time-dependent  calcula¬ 
tions  to  the  mixed  layer  with  parameterization  of  the  interactions 
between  this  layer  and  both  the  underlying  and  overlying  layers. 

Although  the  above  hypotheses  and  assumptions  may  be  con¬ 
sidered  unduly  restrictive,  earlier  work  by  Lavoie  (1972,  197*0 
indicates  that  the  model  provides  a  useful  first  approximation  to 
mesoscale  bounda ry- 1 ayer  weather  patterns  including  convective 
precipitation  induced  by  surface  forcing  factors.  This  justifies 
further  applications  of  the  single  layer  concept  to  those  atmospheric 
situations  where  the  original  hypotheses  and  assumptions  apply. 

Improved  parameterization  of  vertical  fluxes,  better 
vertical  resolution  by  the  introduction  of  more  or  other  variables 
in  the  precipitation  prediction  technique  are  areas  where  changes 
may  be  needed. 


CHAPTER  3 


METHOD  OF  SOLUTION 

3 • 1  Numerical  Technique 

Prognostic  equations  for  0,  h,  and  q  form  a  complete 
system  of  equations  which  may  be  solved  numerically  as  an  initial 
value  problem  with  appropriate  boundary  conditions.  The  forward- 
upstream  finite  difference  technique  described  by  Lavoie  (1972), 
and  used  earlier  by  Estoque  (1962),  Ogura  (1962,  1963)  and  Orville 
(1964,  1965)  among  others,  was  chosen  to  represent  the  partial 
derivatives.  The  scheme  consists  of  (l)  a  forward  time  step, 

(2)  one-sided  spatial  differences  evaluated  upstream  of  each  grid 
point  for  the  advective  terms,  and  (3)  central  differencing  for 
non-advect i ve  terms.  Haltiner  (1971)  has  demonstrated  that  upstream 
differencing  produces  a  damped  solution.  However,  the  damping 
characteristic  of  this  technique  should  not  significantly  modify 
the  results  since  the  disturbances  will  be  continually  reinforced 
and  are  geographically  fixed.  Furthermore,  the  damping  of  short 
wavelengths  should  effectively  prevent  the  nonlinear  numerical 
instibility  that  arises  from  aliasing  of  short  waves.  The 
ac  cy  of  this  scheme  for  simulating  propagating  waves  in 
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comparison  to  other  finite  differencing  techniques  has  been  discussed 
by  Molenkamp  ( 1 968) .  He  indicates  that  it  is  less  accurate  than  some 
alternate  techniques;  however, i t  has  the  advantage  of  requiring  less 
computer  time. 

A  47.6  km  by  47-6  km  grid  mesh  of  1 369  points  was  used  to 
represent  numerically  the  variables  in  the  mesoscale.  The  choice  of 
a  time  increment  was  somewhat  empirical,  because  a  linear  computational 
stability  analysis  of  the  complete  difference  equations  is  very  diffi¬ 
cult.  The  i ncrement  At  was  estimated  simply  from  the  stability  criterion 
for  the  basic  two  dimensional  advection  equation.  Therefore,  the 
stability  criterion  is 

cAt/Ax<2  ^  (Haltiner,  1972) 

where  c  is  the  phase  speed  of  the  fastest  waves,  At  is  the  time 
increment,  and  Ax  is  the  horizontal  space  increment.  A  time  increment 
of  5  min.  appeared  to  be  acceptable  because  no  evidence  of  instability 
showed  up  throughout  the  integration  period. 

Small-scale  perturbations  created  by  fast-moving  gravity 
waves,  boundary  effects  and  truncation  errors  may,  if  left  unchecked, 
'blow-up'  and  obscure  the  meteorological  phenomena.  The  growth  of 
such  spurious  short  waves  was  effectively  controlled  by  explicit 
smoothing  of  the  predicted  fields  after  every  fourth  time  step.  The 
following  smoothing  operator  was  used  to  damp  short  waves 
(Assel in,  1 966)  ; 

s(i,j)  =  A(i,j)  +  b/i)(A(i-1  J)  +  A(i  +  1,j) 

+  A ( i  ,j-1)  +  A(  i  ,  j  +  1 )  -  4A(i  ,j) ) 
where  S(i,j)  is  the  smoothed  grid  point  value,  A(i,j)  is  the  grid 
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point  value  and  b  is  the  smoothing  parameter  set  equal  to  +0.5.  The 
predicted  fields  are  then  'desmoothed'  by  applying  the  above  operator 
a  second  time  with  b  now  equal  to  -0.5*  This  partially  restores  the 
amplitude  of  wavelengths  which  have  been  damped  by  the  original 
smoothing.  The  response  function  of  the  smoothing  and  desmoothing 
operators  is 

1  n 

R(k.)  =  (— )  E  (1  -  b  +  b  cos  k.Ax) 

•  n  .  t  i 

i  =  l 

where  Ax  =  grid  spacing,  b  =  the  smooth  i  ng  parameter ,  n  =  the  number  of 
dimensions  and  k.  is  the  wave  number. 

i 

3 . 2  Numerical  Procedure 

The  computational  steps  required  for  this  numerical  simula¬ 
tion  are  given  below. 

3.2.1  Initial  Condi t i ons 

Wind,  layer  depth,  potential  temperature  and  specific 
humidity  in  the  mixed  layer  are  set  uniform  over  the  entire  grid.  The 
pressure  gradient  is  calculated  assuming  level  terrain  without  heat 
and  moisture  sources.  Friction  typical  of  a  level  land-water  surface 
is  included  so  that  the  wind  vectors  cross  the  isobars  at  a  specified 
angle.  The  Coriolis  parameter  f  is  allowed  to  vary  across  the  region. 

The  procedure  of  calculating  the  pressure  from  the  wind 
field  was  followed  because  the  numerical  solution  of  the  primitive 
equations  is  extremely  sensitive  to  initial  conditions.  The  primitive 
momentum  equation  equates  acceleration  to  the  resultant  of  Coriolis, 


28 


friction  and  pressure-gradient  forces  per  unit  mass.  This  resultant 
is  normally  about  ten  percent  of  the  individual  forces,  so  that  observa¬ 
tional  errors  or  'noise1  due  to  microscale  phenomena  in  the  pressure  or 
wind  fields  may  be  as  large  as  the  accelerations  themselves.  Such 
spurious  accelerations  create  fictitious  gravity  waves  whose  ampli¬ 
tudes  may  be  larger  than  those  encountered  in  the  real  atmosphere 
(Thompson,  1961).  This  problem  can  be  avoided  by  imposing  the  some¬ 
what  artificial  condition  of  balance  between  the  Coriolis,  friction 
and  the  pressure-gradient  forces  at  the  initial  time.  In  view  of  the 
fact  that  the  wind  may  be  considered  the  defining  parameter  for  spring 
upslope  weather,  the  pressure  field  is  calculated  from  the  wind  field; 
this  ensures  a  pressure  field  compatible  with  the  initial  wind  field. 

High  frequency  gravity  waves  in  the  initial  wind  and 
pressure  fields  are  damped  by  applying  the  previously-defined  smooth¬ 
ing  and  desmoothing  operators. 

3.2.2  Readjustment  Process 

At  the  start  of  integration,  the  surface  is  uncovered  and 
the  terrain  relief  is  exposed.  The  momentum  and  heat  fluxes  proceed 
to  alter  the  flow  while,  at  the  same  time,  vertical  motions  induced 
by  the  winds  blowing  against  the  sloping  terrain  begin  to  perturb 
the  upper  boundary  of  the  mixed  layer.  The  readjustment  process 
proceeds  as  follows:  (1)  advance  the  horizontal  wind  components  u 
and  v  to  t  =  At/2,  (2)  use  these  forecast  va  1  ues  of  u  and  v  midway 

through  the  first  time  interval  to  advance  0,  q,  and  h  to  t  =  At, 

(3)  subsequent  u  and  v  steps  are  taken  over  At,  but  the  wind  con¬ 
tinues  to  be  evaluated  at  times  midway  through  the  interval  over 
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which  the  other  three  fields  are  advanced.  This  procedure  permits 

the  maximum  time  step  to  be  nearly  double  that  required  for  linear 

computational  stability  with  the  non-staggered  format  (Lavoie,  1972). 

Integration  is  continued  for  six  hours  of  model  time  to  provide  a 

sufficient  interval  for  the  \^,  h,  0  and  q  fields  to  respond  to  the 

various  forcing  factors.  It  was  found  that  after  six  hours  the  rate 

0 

of  change  of  the  dependent  variables  had  become  sufficiently  small 
that  a  near  steady-state  solution  may  be  assumed. 

3.2.3  Boundary  Conditions 

The  dependent  variables  are  held  constant  at  the  inflow  and 
side  boundaries.  At  the  outflow  boundaries  the  value  of  the  dependent 
variables  are  set  equal,  at  each  time  step,  to  the  values  at  the  next 
grid  point  upstream. 

3 . 3  Grid  Domain,  Topography  and  Drag  Coefficients 

3.3.1  Grid  Domain 

Figure  5  shows  the  area  of  interest  and  the  grid  domain. 

The  grid,  based  on  the  standard  381  km  northern  hemispheric  grid, 
consists  of  37  rows  and  37  columns.  The  spacing  is  equal  to  47.6  km 
in  the  x  and  y  directions,  true  at  60N  latitude  on  a  polar  stereo¬ 
graphic  projection.  The  orientation  of  the  grid  is  such  that  the 
x-axis  is  nearly  parallel  to  the  gradient  of  the  terrain,  with  the 
Canadian  Shield  lying  upstream  of  the  Western  Plains  for  an  east- 
northeasterly  flow.  The  Continental  Divide  provides  a  natural 
boundary  to  the  southwest,  while  Hudson  Bay  serves  a  similar  function 
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Fig.  5. 


Grid 


doma i n 


(The  outermost  ring  of  intersections 
represents  the  grid  boundary.) 
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to  the  northeast.  The  northwest  and  southeast  edges  of  the  grid  were 
chosen  to  incorporate  the  maximum  amount  of  the  Canadian  Western  Plains. 

3.3.2  Topography 

The  terrain  elevation  at  each  grid  point  was  obtained  from 
topographic  maps^  of  scale  1:250,000.  A  large  portion  of  these  data  had 
been  abstracted  previously  by  Charette  (1971).  The  terrain  grid-point 
values  were  smoothed  using  the  same  smoothing  operator  that  was  sub¬ 
sequently  applied  to  the  predicted  fields.  This  provided  a  smoothed 
topographical  field  on  the  same  scale  as  the  atmospheric  motion  being 
cons i dered . 

3.3.3  Drag  Coefficients 

In  order  to  include  fluxes  of  momentum,  sensible  heat  and 
moisture  in  the  model,  it  is  necessary  to  know  the  spacial  distribution 
of  Cp,  the  drag  coefficient,  C^,  the  sensible  heat  transfer  coefficient 
and  C^,  the  evaporation  coefficient  over  the  area  of  interest.  It  is 
also  desirable  to  relate  these  coefficients  to  the  horizontal  wind 
speed  and  to  atmospheric  stability. 

Following  Cressman  (i960),  we  discriminate  between  the  drag 
coefficient  over  flat  land  or  water  (skin  drag),  and  the  drag  due  to 
terrain  relief  (terrain  drag).  The  drag  coefficient  CD  is  considered 
to  be  made  up  of  two  partial  drag  coefficients  C^  and  C^,  where  C^  is 


Topographic  series  produced  and  printed  in  Canada  by  the 
Surveys  and  Mapping  Branch,  Department  of  Energy,  Mines  and  Resources, 
Government  of  Canada  and  in  the  United  States  by  the  Interior  Geo¬ 
logical  Survey,  Washington,  D.C. 


32 


the  skin  drag  and  is  the  terrain  drag  coefficient.  Values  for 
for  the  standard  38 1  km  northern  hemispheric  grid  were  supplied  by 
Dr.  Cressman  of  NOAA  through  G.  R.  Schram.  These  values  were  extended 
to  the  A7.6  km  grid,  a  subset  of  the  standard  grid,  by  linear  inter¬ 
polation  with  minor  subjective  changes.  Rather  than  using  a  constant 
value  of  C1  as  Cressman  (I960)  has  done,  an  attempt  was  made  to  relate 
the  skin  drag  coefficient  to  the  horizontal  wind  speed  (m  sec"1)  using 
the  following,  formula  (Deacon  and  Webb,  1 962 ) : 


Z  =  (1  +  0 . 7 | V | ) 0. 001 

Given  the  drag  coefficient  for  the  neutral  atmosphere  C^n 
as  a  function  of  wind  speed  and  terrain  relief,  we  proceed  to  adjust 
these  values  according  to  atmospheric  stability  after  Deardorff  (1968). 
For  the  unstable  case,  which  is  the  only  one  that  applies  in  this 
simulation,  the  drag  coefficient  becomes: 


C0  =  (1  +  (7/a)  Log  (1  -  aR.b))CDn 


n  ? 

where  a  =  0.83  "  '  ,  and  R.,  is  the  bulk  Richardson  number,  a 

Dn  ib 

measure  of  atmospheric  stability. 
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where  Z  is  the  height  where  V,  0  and  q  are  known,  g  is  the  acceleration 

due  to  gravity,  the  subscript  o  refers  to  surface  values  with  | Vq |  =  0 

and  T  is  the  virtual  temperature,  here  assumed  to  be  equal  to  T  which 
v 

is  given  by  (T^  +  Tq)/2. 

The  sensible  heat  and  water  vapour  transfer  coefficients 
are  obtained  by  assuming  CH  =  =  CD  f°r  the  neutra1  case  and  then 


adjusting  the  values  of  C  on  the  basis  of  stability  after  Deardorff 

H 

(1968),  while  continuing  to  assume  that  C..  =  C  .  Therefore,  the  heat 

t  H 

and  water  vapour  transfer  coefficients  become: 

CH  =  CE  =  0  +  (11/c)  Log  (1  -  cRjb))  C0n 

where  c  =  0.25  Cn  0.80. 

Dn 

The  fields  of  drag  coefficient  and  transfer  coefficient  grid- 
point  values  are  smoothed  using  the  previously  defined  smoothing  and 
desmoothing  operators.  These  values  are  assumed  to  remain  constant 
during  the  integrating  period. 
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CHAPTER  A 


EXPERIMENTAL  DESIGN,  INITIALIZATION,  AND  RESULTS 

4 .  1  Introductory  Comments 

It  has  been  hypothesized  after  Lavoie  (1972,  197*0  that  a 
single  well-mixed  layer  may  be  used  to  represent  many  of  the  signifi¬ 
cant  aspects  of  mesoscale  bounda ry- 1 ayer  weather  disturbances  and 
that  these  disturbances  provide  control  over  the  location  and  intensity 
of  convective  activity.  The  present  model  has  been  designed  to  study 
the  evolution  of  a  well-mixed  planetary  boundary  layer  in  response  to 
several  forcing  factors  including  terrain  rel i ef , . f r i ct ion ,  surface 
heating  and  the  release  of  latent  heat  accompanying  the  formation  of 
precipitation.  This  represents  an  initial  attempt  to  simulate 
quantitatively  the  temporal  and  spatial  evolution  of  mesoscale  atmos¬ 
pheric  weather  patterns  resulting  from  the  interaction  of  an  east- 
northeasterly  circulation  of  moderately  moist  arctic  air  with  the 
relatively  warm,  sloping  Canadian  Western  Plains  during  late  spring. 

Vertically  homogeneous  initial  distributions  of  the 
dependent  variables  layer  depth,  wind,  specific  humidity  and  potential 
temperature  were  used.  Detailed  upstream  topography  was  not  delineated 
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initially,  because  there  is  little  relief  over  the  Canadian  Shield. 
Friction  was  assumed  to  be  the  same  over  the  whole  grid  so  as  to 
provide  balanced  flow  with  respect  to  the  Shield  at  all  grid  points. 

At  t  =  0  the  Western  Plains  was  'uncovered'  with  its  topography, 
roughness  differential,  moisture  sources  and  stronger  surface  heating. 
The  readjustment  procedure  was  followed  at  each  grid  point  by  means 
of  equations  (19),  (23),  (30),  and  (3^0  using  the  forward  upstream 
numerical  integration  scheme.  Integration  was  continued  for  six  hours 
of  model  time  to  provide  a  sufficient  interval  for  the  PBL  to  respond 
to  the  various  surface  forcing  factors.  It  was  found  that  after  six 
hours  the  rate  of  change  of  the  dependent  variables  had  become  suffic¬ 
iently  small  that  a  near  steady  state  solution  was  achieved. 

Lavoie  (1972,  197^0  has  shown  that  the  deformation  of  the 
inversion  surface  from  its  initial  height  is  a  useful  guide  to  the 
intensity  and  distribution  of  mesoscale  convective  precipitation. 
Following  Lavoie  but  on  a  scale  nearly  one  order  of  magnitude  larger 
than  his,  we  proceed  using  the  hypothesis  that  the  atmospheric  re¬ 
sponse  on  the  mesoscale  is  dictated  by  the  behaviour  of  the  mixed 
layer  as  a  whole  with  only  secondary  dependence  on  the  detailed  vertical 
structure. 

k . 2  Simulation  of  Typical  Late  Spring  Upslope  Weather 

The  most  representative  low-level  flow  conditions  for  late 
spring  over  the  Canadian  Western  Plains  are  similar  to  the  mean,  for 
during  the  encompassing  period,  March  through  May,  the  average  surface 
circulation  persists  out  of  the  east-northeast  for  all  but  the  western 
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regions  (Bryson,  1966).  Typical  daytime  values  of  atmospheric 
structure  and  low  level  flow  were  selected  to  test  the  ability  of  the 
model  to  simulate  late  spring  upslope  weather.  The  term  upslope  is 
usually  applied  to  a  shallow  east-northeasterly  circulation  associated 
with  outbreaks  of  arctic  air  masses  over  the  Western  Plains.  Such 
easterlies  are  usually  confined  to  the  lowest  1  to  2  km  of  the 
atmosphere. 

During  the  spring  thaw  period  the  surface  of  the  Canadian 

Shield,  which  is  normally  over  50  percent  water  in  the  summer-time, 

probably  approaches  a  true  water  or  'lake'  surface  with  a  temperature 

near  0C.  Continental  Arctic  air  (cA)  is  modified  rapidly  to  a  brand 

of  air  very  similar  in  character  to  a  winter  cold  maritime  Arctic 

(cmA)  air  mass,  with  typical  wet-bulb  potential  temperatures  in  the 

range  -5C  <0  <  -1C.  This  modified  arctic  air  is  usually  referred 

w 

to  as  transitional  continental  Arctic  (tcA) .  The  resulting  unstable 
turbulent  low-level  flow,  further  perturbed  by  upslope  and  stronger 
surface  heating  upon  reaching  the  Western  Plains  often  produces 
extensive  decks  of  limited  convective  clouds  (stratocumul us) .  Pre¬ 
cipitation  will  generally  be  light  and  intermittent  and  may  take  the 

form  of  drizzle  or  snow.  The  duration  of  these  conditions  may  be 

/ 

several  days. 

At  time  t  =  0  an  east-northeasterly  flow  of  arctic  air, 
modified  to  tcA  by  interacting  with  the  Shield,  was  assumed  to  cover 
the  entire  grid  domain.  The  model  air  mass  was  horizontally  uniform 
with  the  same  vertical  structure  and  low-level  flow  at  all  grid  points. 
For  t  >  0  the  topography  was  uncovered.  As  the  integration  proceeded 
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mesoscale  weather  patterns  with  areas  of  limited  convective  precipita¬ 
tion  evolved  due  to  spatial  and  temporal  changes  in  the  well-mixed 
boundary  layer. 

The  'undisturbed'  air  mass  and  surface  parameters  used  to 

initialize  the  simulation  of  'typical'  late  spring  upslope  weather 

are  outlined  below.  Fig.  6  depicts  the  tcA  sounding.  The 

wet  bulb  potential  temperature  0  was  taken  to  be  -3C,  the 

w 

surface  temperature  Tq  was  assumed  to  be  OC ,  and  the  inversion 
height  h  was  set  at  1.4  km,  all  in  accordance  with  the  previous 
discussion.  The  thermal  structure  of  the  model,  the  Normand  principle 
(Normand,  1938),  and  the  above  values  determined  the  surface  dew¬ 
point  to  be  near  -6.7C.  The  surface  height  of  the  Shield,  Zq,  was 
set  at  250  m  above  ms  1 ;  detailed  terrain  was  ignored  at  the  initial 
time  because  there  is  little  relief  over  the  Shield.  The  vertical 
structure  or  representative  sounding  of  the  model  atmosphere  was  then 
built  up  from  surface  values.  This  was  done  because  a  predictive  role 
using  actual  data  was  envisaged  in  the  model's  future  and  it  was  felt 
desirable  to  build  into  the  model  the  ability  to  use  the  wealth  of 
surface  data  rather  than  the  relatively  scarce  upper  air  data.  The 
lapse  rate  of  temperature  was  arbitrarily  set  at  10. 2C  km  1  for  the 

constant-flux  layer  from  the  surface  Z  to  the  anemometer  level  Z  , 

o  s 

and  a  dry-adiabatic  lapse  rate  of  9 . 8C  km  1  was  required  by  the  model 
atmosphere  from  Z^  to  h,  the  top  of  the  well-mixed  layer.  A  near 
isothermal  layer  was  assumed  above  h.  The  superad i abat i c  surface 
layer  ensured  an  initial  upward  heat  flux.  The  lapse  rate  of  dew¬ 
point  was  set  at  10. 1C  km  1  from  the  surface  to  level  h.  This 
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decrease  which  is  slightly  in  excess  of  the  decrease  in  temperature 

with  height  maintained  a  nearly  constant  0  value  through  the  boundary 

w 

layer  and  it  also  ensured  an  initial  upward  flux  of  moisture.  The 
above  values  were  then  used  to  compute  the  model  parameters  of  potential 
temperature  and  specific  humidity  at  the  various  levels  (see  Appendix). 
The  cloud  base  given  by  the  lifting  condensation  level  (LCL)  ,  or  the 
level  to  which  unsaturated  air  would  have  to  be  raised  in  a  dry- 
adiabatic  expansion  to  produce  condensation,  was  assumed  uniform 
throughout  the  air  mass  and  was  calculated  as  outlined  in  Pettersen 
(1956)  (see  Appendix). 

The  greater  proportion  of  land  to  water  surface  on  the 
Western  Plains  coupled  with  the  earlier  loss  of  snow  cover  allows 
the  surface  of  the  plains  to  heat  under  day-time  conditions  to  a 
potential  temperature  in  excess  of  the  surface  potential  temperature 
0q  on  the  Shield.  An  excess  surface  potential  temperature  0q+  value 
of  kK  was  chosen  to  represent  mid-day  conditions  on  the  plains.  A 
transition  area  between  the  Shield  and  the  plains  was  included  with 
surface  excess  potential  temperature  taken  to  be  0q+/2. 

Since  vertical  homogeneity  of  the  horizontal  wind  was 
assumed,  a  representative  mixed-layer  wind  was  taken  from  the  mean 
surface  (anemometer  level)  streamline  pattern  for  May  given  by 
Bryson  (1966),  Fig.  1.  This  is  the  mean  flow.  The  periods  of 
interest  are  those  days  when  the  discontinuity  in  the  wind  field 
lies  just  east  of  the  Rockies  subjecting  the  entire  Western  Plains 
area  to  an  east-northeasterly  flow.  A  horizontal  wind  from  075 
degrees  at  7  rn  sec  1  was  taken  as  representative  of  conditions  over 


the  Shield.  This  was  assumed  to  be  the  flow  at  all  grid  points 
initially.  A  detailed  wind  field  was  not  included  initially  for  two 
reasons:  (a)  it  is  impossible  to  obtain  a  horizontal  wind  field 

that  has  not  already  been  influenced  by  topography  and  (b)  using  a 
uniform  field  avoided  the  problem  of  including  real  or  fictitious 
lines  of  convergence  or  divergence  at  the  initial  time  which  might 
mask  the  'true'  evolution  of  the  horizontal  wind  field  in  response 
to  the  various  surface  forcing  factors. 

The  horizontal  pressure  gradients  were  calculated  from  the 
wind  field  to  avoid  fictitious  acce 1 erat i ons  as  discussed  in  Chapter  3. 
The  details  of  these  calculations  are  outlined  in  the  Appendix. 

4.2.1  Resul ts 

Time-mean  precipitation  maps  based  on  the  sub-synoptic  net¬ 
work  of  climatic  stations  show  the  average  distribution  of  precipita¬ 
tion  in  the  mesoscale.  Cyclonic  and  purely  convective  precipitation 
are  expected  to  be  comparatively  irregular  in  their  incidence  and 
their  effects  become  smoothed  when  long  periods  are  considered.  On 
the  other  hand,  surface  forcing  factors  or  topographical  influences 
are  permanent  features.  They  produce  local  deviations  from  the 
regional  average  precipitation  and  their  effects  stand  out  on  most 
monthly  or  longer-period  precipitation  maps  (Douglas  and  Glaspoole, 
1947) •  Thus,  a  good  estimate  of  the  topographical  precipitation  of 
a  region  may  be  obtained  from  long-term  climati  records  (Reinelt, 
1969).  The  mean  (1931"60)  precipitation  map  for  the  encompassing 
spring  period,  March  to  May,  clearly  shows  that  the  higher  elevations 
in  the  Western  Plains  also  have  higher  precipitation  (Longley,  1972). 
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The  mean  spring  precipitation  patterns,  therefore,  were  chosen  to 
verify  the  present  simulation  of  late  spring  upslope  weather  (see 
Fig.  7).  The  period  which  has  been  referred  to  as  late  spring  varies 
somewhat  from  year  to  year  with  respect  to  the  calendar,  thus  the 
March  to  May  precipitation  map  was  considered  most  appropriate  for 
verification  purposes. 

The  topographic  precipitation  component,  i.e.,  precipitation 
due  to  surface  forcing  factors,  for  the  Western  Plains  can  be  estima¬ 
ted  if  some  reference  precipitation  value  can  be  determined  which  is 
devoid  of  all  topographic  influences.  Figure  8  shows  the  precipitation 
component  due  to  surface  forcing  factors  based  on  the  hypothesis  that 
the  generally  observed  3  inch  value  of  precipitation  is  entirely  due 
to  non- topog raph i c  effects.  This  leaves  10-30  percent  of  the  total 
precipitation  attributable  to  surface  forcing  factors  in  agreement 
with  Reinelt  (1969)  that  topographic  precipitation  represents  about 
37  percent  of  the  total  for  the  foothills  and  more  than  13  percent 
for  the  rest  of  the  Western  Plains.  Figures  7  and  8  show  precipita¬ 
tion  maxima,  or  excesses  due  to  surface  influences,  in  two  northwest- 
southeast  bands.  One  band  lies  over  the  foothills  with  fingers 
extending  over  the  Swan  and  Cypress  hills  and  a  second  band  lies 
along  the  Manitoba  escarpment  which  includes  the  Turtle,  Riding  and 
Duck  mountains,  as  well  as  the  Pembina,  Porcupine,  and  Pasquia  hills. 
This  second  band  extends  northwestward  from  the  escarpment  just  south¬ 
west  of,  and  paralleling,  the  edge  of  the  Canadian  Shield.  The  north¬ 
west-southeast  band  structure  of  the  precipitation  pattern,  as  well 
as  the  dry  pocket  immediately  west  of  the  Duck  Mountains  and  Riding 
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Mean  March  to  May  precipitation  (inches)  (after  Longley,  1972) 
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Mountains,  clearly  shows  the  effect  of  the  east-northeasterly  circula¬ 
tion  of  moderately  moist  arctic  air  ascending  the  upland  areas  of  the 
Western  Plains  during  the  spring  period. 

Figure  9  shows  the  simulated  precipitation  pattern  produced 
by  the  model  using  'typical'  late  spring  upslope  weather  conditions. 
Although  the  values  of  the  simulated  precipitation  rate  and  the  mean 
topographical  precipitation  for  three  months  cannot  readily  be  com¬ 
pared  quantitatively,  it  can  be  seen  that  the  simulated  wet  and  dry 
patterns  correspond  quite  well  with  the  actual  topographic  precipita¬ 
tion  patterns.  A  number  of  factors  are  believed  to  contribute  to 
the  difference  in  patterns.  These  include  (a)  deviations  of  actual 
meteorological  conditions  from  those  used  to  initialize  the  simulation, 
(b)  the  simulation  included  day-time  heating,  whereas  the  mean  pre¬ 
cipitation  pattern  includes  both  day  and  night  precipitation,  (c)  the 
lack  of  mesoscale  detail  available  for  the  drag  and  transport  co¬ 
efficient  distributions,  and  (d)  evaporation  of  precipitation  below 
cloud  base,  and  wind  transport  of  droplets  have  been  neglected.  It 
was  also  realized  that  the  spacing  and  location  of  climatic  stations 
lead  to  inaccuracies  in  the  resultant  mean  precipitation  pattern. 
Nevertheless,  considering  the  above,  the  simulated  precipitation 
pattern  is  remarkably  similar  to  the  mean  spring  topographic  pre- 
ci pi  tat  ion  pattern . 

\ 

Figure  10  shows  the  vertical  velocities  generated  after 
six  hours  of  model  time.  The  large-scale  relief  of  the  Prairie 
Provinces  may  be  visualized  as  a  gently-sloping  wedge  rising  from 
the  Hudson  Bay  lowlands  in  the  northeast  to  the  Rocky  Mountains  in 
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Fig.  9.  Simulated  topographic  precipitation  pattern  (inches  day  1). 
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Simulated  vertical  velocity  pattern  (cm  sec  1). 


the  west-southwest.  An  east-northeaster  1 y  circulation  will,  therefore, 
experience  forced  ascent  and,  as  a  result,  small  positive  vertical 
velocities  are  expected  over  the  entire  region.  Cross  isobaric  flow 
leading  to  horizontal  divergence,  and  the  resultant  synoptic-scale 
subsidence  will  add  a  negative  component  to  the  vertical  velocity 
values  as  the  wind  turns  more  and  more  out  from  the  high  pressure 
centre  in  the  northeast  due  to  increasing  surface  roughness  from 
the  lowlands  to  the  Rockies. 

The  simulated  vertical  velocity  pattern  in  the  mesoscale 
reveals  other,  generally  more  significant,  features.  Looking  from 
east  to  west  in  Figure  10, one  sees  an  area  of  marked  subsidence  where 
the  Lake  of  the  Woods  area  drops  off  to  Lake  Winnipeg.  Further  west, 
there  is  a  north-south  band  of  relatively  strong  ascent  along  and 
northwest  of  the  Manitoba  escarpment  with  maxima  immediately  upwind 
of  the  individual  upland  areas.  Corresponding  cellular  areas  of 
subsidence  or  reduced  ascent  are  observed  immediately  downwind  of 
these  upland  areas.  The  extreme  vertical  velocities  in  these  cells 
are  listed  in  Table  1.  The  general  area  of  positive  vertical  velociti 
associated  with  the  escarpment  can  be  attributed  to  the  abrupt  upslope 
of  the  relief.  The  cellular  maxima  and  minima  of  vertical  velocities 
associated  with  the  individual  upland  areas  are  due  to  the  slowing  of 
the  low  level  flow  as  it  ascends  the  ridges  and  a  corresponding  in¬ 
crease  in  wind  speed  immediately  downstream  from  these  areas.  With 
the  initial  wind  direction  nearly  perpendicular  to  the  orientation 
of  the  escarpment,  significant  areas  of  resultant  convergence  upstream 
and  divergence  downstream  from  the  individual  hills  produced  the 
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Table  1 


Extreme  vertical  velocities  -  Manitoba  Escarpment 


Upland  Area 

Max.  Upwind 

Min.  Downwind  W, 

h 

Riding  Mountains 

+2.1 

cm 

sec*  1 

-0.9 

cm 

sec"  1 

Duck  Mountains 

+2.7 

cm 

sec  1 

-1.3 

cm 

sec  1 

Porcup i ne  Hills 

+2.3 

cm 

sec"  1 

-0.8 

cm 

sec  1 

Pasqu i a  Hills 

+2.1 

cm 

sec"  1 

-0.4 

cm 

sec*  1 

Pemb i na  Hills 

+2.4 

cm 

sec  1 

-0.9 

cm 

sec”  1 

Wapawekka  Hills 

+2.9 

cm 

sec  1 

-0.5 

cm 

sec"  1 

observed  maxima  and  minima  in  the  vertical  velocity  pattern. 

The  plains  region  contains  a  number  of  other  significant 
upland  areas  each  of  which  acts  through  forced  ascent  and  frictional 
effects  to  produce  vertical  velocities  as  shown  in  Table  2.  The 
foothills  in  the  west  produced  a  north-south  band  of  positive  vertical 
velocities  with  a  maximum  value  of  2.4  cm  sec  *. 

The  vertical  velocity  patterns  described  above  are  further 
complicated  by  topographic  deflection  which  acts  to  produce  local 
areas  of  convergence  on  the  upwind  shoulders  of  the  hills  and  areas 
of  divergence  on  the  downwind  shoulders. 

Thus,  it  can  be  seen  that  the  effects  of  topographical 
influences  on  the  mesoscale  vertical  velocity  pattern  are  far  more 
intricate  than  just  general  upslope  with  the  resulting  vertical 
velocity  values  reduced  by  synoptic  scale  subsidence.  The  general 


Table  2 


Extreme  vertical  velocities  -  other  upland  areas 


Upland  area  Max.  Upwind  Min.  Downwind 


Moose  Mountains 

+2.0 

cm 

sec  1 

-0.6 

cm 

sec”1 

Touchwood  Hills 

+  1.5 

cm 

sec  1 

-0.5 

cm 

sec  1 

Dirt  Hills 

+  1.5 

cm 

sec  1 

-0.5 

cm 

sec  1 

Cypress  Hills 

+  1.7 

cm 

sec  1 

-0.4 

cm 

sec  1 

Cheecham  Hills 

+  1.2 

cm 

sec  1 

-0.7 

cm 

sec  1 

Swan  Hills 

+2.8 

cm 

sec  1 

-1.2 

cm 

sec  1 

B i rch  Mounta i ns 

+3.5 

cm 

sec”1 

-2.3 

cm 

sec” 1 

upslope  motion  is  complicated  by  mesoscale  frictional  convergence 
and  divergence  effects  associated  with  the  individual  upland  areas, 
subsidence  resulting  from  local  downslope  and  the. effect  of  spatial 
variations  of  momentum  fluxes  due  to  nonuniform  surface  roughness. 

Figure  11  shows  the  corresponding  inversion  deformation 
pattern.  It  has  a  marked  similarity  to  the  terrain  relief  with  the 
displacement  of  the  top  of  the  boundary  layer  or  inversion  base  havi 
maxima  generally  one  grid  length  upwind  from  the  ridges  or  upland 
peaks.  Defining  the  location  of  an  upland  area  as  the  position  of 
the  appropriate  local  maximum  grid  point  terrain  height  value,  the 
inversion  deformation  was  found  to  vary  from  20  to  50  percent  of  the 
terrain  height  as  shown  in  Table  3* 
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Simulated  inversion  deformation  pattern  (m) . 
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Table  3 

Inversion  deformations 


Upland  area 

Ter ra i n 

Ht.  (m) 

1 nvers i on 
Deformation  (m) 

Ah/Z  x  100 
o 

Z 

Ah 

o 

Footh ills* 

983 

399 

40.6 

Cypress  Hills 

1086 

375 

34.5 

Swan  Hills 

1 084 

434 

40.0 

Dirt  Hills 

926 

276 

29.8 

Moose  Mountains 

723 

249 

34.4 

Touchwood  Hills 

684 

248 

36.3 

B i rch  Mounta i ns 

702 

259 

36.9 

Cheecham  Hills 

691 

151 

21.9 

Riding  Mountains 

632 

245 

38.8 

Duck  Mountains 

639 

304 

47.6 

Porcupi ne  Hills 

596 

276 

46.3 

Pasqu i a  Hills 

523 

237 

45.3 

Waspaweki a  Hills 

597 

204 

34.2 

Turt 1 e-Pemb i na 

589 

271 

46.0 

*The  foothills  are  defined  as  a  test  area  of  14  grid 
points  west  of  Red  Deer. 
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Due  to  the  size  and  detail  of  the  area  of  interest,  a 
selected  cross  section,,  approximately  parallel  to  the  flow,  was  con¬ 
structed  to  display  graphically  the  relationship  between  the  terrain 
relief,  vertical  velocities,  inversion  deformation,  and  limited 
convective  precipitation.  This  cross  section,  shown  in  Figs.  12,  13 
and  14,  reveals  an  apparent  inconsistency.  The  vertical  velocity 
associated  with  the  Riding  Mountains  was  larger  than  that  associated 
with  the  foothills,  yet  the  inversion  deformation  and  the  rate  of 
precipitation  was  greater  in  the  latter  area.  This  resulted  from  the 
much  larger  transfer  coefficients  and  consequently  greater  vertical 
heat  fluxes  in  the  foothills  region.  The  increase  in  the  thickness 
of  the  mixed  layer  was  due  to  the  last  term  of  equation  (30).  The 
increase  resulted  from  the  incorporation  of  air,  whose  potential 
temperature  was  the  same  as  that  of  the  column,  in  order  to  prevent 
the  development  of  a  superad iabat i c  layer  at  the  top  of  the  mixed 
layer  whenever  the  capping  inversion  was  wiped  out  by  heating  from 
below.  This  simulates  the  well-known  effect  that  mountains  and  their 
foothills  act  as  effective  heat  sources  which  tend  to  promote  precipita¬ 
tion  of  a  convective  nature.  The  subsequent  release  of  latent  heat 
further  warms  the  layer  and  increases  the  inversion  deformation.  Thus, 
it  can  be  seen  that  topography  contributes  significantly  both  to 
thermal  convection  and  to  forced  ascent  with  frictional  effects.  How¬ 
ever,  due  to  mutual  interactions  it  is  usually  difficult  to  separate 
the  relative  contributions  of  the  various  surface  forcing  factors  to 


the  resultant  limited  convective  showers. 


Height  (km) 
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Vertical  velocity  -  topography  cross-section 
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Inversion  deformation  -  topography  cross-section 
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Fig.  1 b . 


Precipitation  -  topography  cross-section. 
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^4 . 3  A  Sensitivity  Experiment 

The  relative  contribution  of  surface  heating  compared  with 
the  topographical  effects  of  upslope  and  frictional  convergence  or 
divergence  to  the  evolution  of  mesoscale  bounda ry- 1 ayer  weather 
patterns  with  limited  convective  precipitation  was  investigated  by 
running  the  model  wh i 1 e-  neg 1 ect i ng  surface  eddy  heat  fluxes.  Meteoro¬ 
logical  conditions  identical  to  those  used  in  the  preceding  simulation 
were  used  so  that  the  various  resultant  fields  could  be  compared. 

The  sensitivity  run  was  an  attempt  to  isolate  the  influence  of  surface 
heating  and  to  assign  a  relative  magnitude  to  it.  Although  the  various 
surface  forcing  factors  are  highly  i nterrel ated, some  separation  was 
possi ble. 

The  experiment  involved  setting  the  transfer  coefficients 
for  heat  equal  to  zero  throughout  the  integration  period.  Vertical 
fluxes  of  momentum  and  moisture  were  not  interrupted. 

Figure  15  shows  the  inversion  deformation  generated  after 
six  hours  of  model  time.  The  most  obvious  differences  between  this 
pattern  and  the  inversion  deformation  pattern  obtained  in  the  simula¬ 
tion  of  typical  late  spring  upslope  weather  were  the  generally  much 
smaller  values  of  the  inversion  deformation.  The  relative  magnitude 
of  the  reduction  in  the  inversion  deformation  may  be  illustrated 
by  considering  the  two  equal-size  test  areas  outlined  in  Figure  13. 

The  removal  of  vertical  eddy  heat  fluxes  resulted  approximately  in  a 
35  percent  reduction  in  the  inversion  deformation  over  the  western 
area  where  the  transfer  coefficients  were  relatively  large  and  a  10 
percent  reduction  for  the  eastern  area  where  the  transfer  coefficients 
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Simulated  inversion  deformation  pattern, 
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were  smaller  in  magnitude.  This  indicates  the  relative  importance  of 
surface  heating  to  the  resulting  inversion  deformation  and  ultimately 
to  the  intensity  of  the  limited  convective  precipitation.  Figure  16 
shows  the  reduced  precipitation  area  that  resulted  when  vertical  eddy 
heat  fluxes  were  neglected.  This  further  illustrates  the  significance 
of  surface  heating.  Thus,  it  was  concluded  that  the  area  and  intensity 
of  limited  convective  precipitation  which  occurs  when  a  moderately 
moist  unstable  arctic  air  mass  advances  upslope  over  the  Western  Plains 
is  significantly  enhanced  by  the  inclusion  of  surface  heating. 

b . b  Simulation  of  an  Observed  PBL  Weather  Pattern 

The  model  has  been  developed  through  consideration  of  con¬ 
ditions  characteristic  of  late  spring  upslope  weather  on  the  Canadian 
Western  Plains.  The  following  case  study  provided  an  independent 
test  of  the  ability  of  the  model  to  simulate  boundary-layer  weather 
patterns  derived  from  observations. 

In  order  to  reproduce  a  persistent  limited  convective 
precipitation  situation  maintained  by  an  east-northeasterly  circulation 
of  tcA  air  interacting  with  the  Western  Plains,  we  began  by  establish¬ 
ing  representative  air  mass  and  flow  parameters  for  the  period  and 
area  of  interest.  This  was  achieved  through  an  examination  of  the 
tephigrams  and  hodographs  for  the  upper  air  stations  within  or  near 
the  grid.  Representative  soundings  were  then  chosen  and  used  to 
establish  the  air  mass  structure  specified  to  be  uniform  over  the 
entire  grid  at  t  =  0.  A  low-level  flow  pattern  was  obtained  by 
calculating  the  vertical  mean  of  the  mixed-layer  winds  for  the 
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Fig.  16.  Simulated  precipitation  pattern, 
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representative  upper  air  stations. 

4.4.1  Case  Study  (April  26,  1969) 

The  synoptic  situation  of  April  26,  1200Z  (Figs. 17,  1 8 ,  19) 
has  been  described  by  Reinelt  (1969)  as  being  representative  of 
shallow  upslope  circulations  east  of  the  Rocky  Mountains.  The  surface 
map  has  an  Arctic  high  pressure  centre  northeast  of  Great  Bear  Lake 
with  a  ridge  line  extending  southeastward  over  Hudson  Bay.  The  high, 
in  conjunction  with  a  deep  low  located  in  North  Dakota,  has  established 
a  persistent  northeasterly  flow  over  most  of  the  Western  Plains.  A 
surface  trough,  which  extended  northwestward  from  the  low-pressure 
centre  along  a  line  from  the  Cypress  Hills  to  just  east  of  the  Swan 
Hills,  has  produced  a  northwesterly  flow  over  Alberta  west  of  the 
trough  line.  The  low-level  flow  was  remarkably  similar  to  the  mean 
May  surface  streamline  pattern  (Fig.  l).  The  700  and  500  mb  maps 
(Figs.  18  and  19)  showed  a  moderate  to  strong  northeasterly  circula¬ 
tion  of  considerable  depth  in  Saskatchewan  and  Manitoba  north  of 
the  intense  low,  but  over  Alberta  the  northeasterly  flow  had  become 
weak  to  moderate  and  quite  shallow  under  the  ridge  that  was  intensify¬ 
ing  at  upper  levels. 

A  narrow  spur  of  precipitation  paralleling  the  Rockies  was 
maintained  by  upslope  flow  (Reinelt,  1969),  and  was  intensified  by 
low-level  convergence  along  the  trough.  A  second  area  of  precipitation 
in  southern  Manitoba  and  Saskatchewan  was,  at  least  partly,  of  cyclonic 
origin.  The  trace  of  precipitation  in  northern  Manitoba  was  attrib¬ 
uted  to  weak  upslope  and  to  low  level  convergence  due  to  onshore 


flow. 
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.  17*  Surface  analysis  261200  GMT 
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Fig.  18. 


700  mb  analysis  261200  GMT,  April 


1969. 
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Fig.  19.  500  mb  analysis  261200  GMT,  April  ’969. 
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An  examination  of  the  tephigrams  for  the  upper  air  stations 
within  or  near  the  grid  indicated  that  the  261200  GMT  sounding  for 
The  Pas  was  representative  of  conditions  in  the  arctic  air  mass  modified 
to  tcA  by  interacting  with  the  Canadian  Shield,  but  not  yet  modified 
by  contact  with  the  Western  Plains.  This  was  the  desired  state  for 
initialization  of  the  model.  Air  mass  parameters  for  the  model, 
therefore,  were  taken  mainly  from  the  261200  GMT  sounding  for  The  Pas. 
Some  adjustments  were  made  based  on  the  270000  GMT  sounding  for  The 
Pas  and  the  261200  GMT  and  270000  GMT  soundings  for  Fort  Smith, 

Edmonton  and  Glasgow  in  order  to  obtain  the  most  representative  con¬ 
ditions  in  the  tcA  air  mass.  The  initial  low-level  flow  east  of  the 
trough  was  obtained  by  calculating  the  vertical  mean  wind  from  The 
Pas'  261200  GMT  hodograph,  using  only  those  levels  which  were  below 
the  inversion.  A  similar  procedure,  based  mainly  on  Edmonton's 
hodographs,  gave  the  representative  low-level  flow  west  of  the  trough. 
This  provided  the  flow  pattern  over  the  area  of  interest.  The  initial 
conditions  were:  (a)  Tq  =  0C,  (b)  TdQ=-6C,  (c)  h  =  1350  m,  and 

(d)  V  =  050/8  m  sec*1  east  of  the  trough  and  335/8  m  sec  1  west  of 
the  trough,  where  the  trough  was  simulated  by  a  column  three  grid 
points  wide  over  which  a  linear  transition  in  wind  direction  from 
335  to  050  degrees  took  place.  The  snow-line  for  261200  GMT  is 
included  in  Figure  17.  The  areas  covered  with  snow  were  assumed 
to  have  surface  temperatures  of  0C  as  was  the  Shield  in  the  east 
and  the  mountains  along  the  western  edge  of  the  grid.  The  remaining 
area  was  assumed  to  have  an  excess  surface  potential  temperature  of 
4K  due  to  a  smaller  albedo.  Narrow  transition  zones  were  included  as 
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in  the  previous  simulation.  As  no  clear-cut  outflow  boundary  existed, 
the  dependent  var i abl es . were  held  constant  at  all  lateral  boundaries. 

The  upstream  and  downstream  soundings  (Figs.  20,  21,  22,  23) 
for  the  period  of  the  weather  disturbance  show  the  persistence  of  the 
air  mass  structure  being  considered  as  well  as  some  of  the  changes 
that  occurred  to  the  inversion  height  over  time  and  space.  The  low- 
level  wind  profiles  (see  Table  4)  for  the  same  stations  indicated 
that  the  assumption  of  vertical  homogeneity  in  the  horizontal  wind 
field  was  reasonable. 

A  successful  simulation  was  defined  as  the  reproduction  of 
the  precipitation  pattern  which  actually  occurred  on  April  26,  1969. 
For  verification  purposes  the  actual  precipitation  pattern  was 
obtained  from  the  daily  records  of  the  climatic  stations  within  the 
area  of  interest  (Fig.  24). 

4.4.2  Resul ts 

Due  to  the  nature  of  this  case,  only  the  spur  of  precipita¬ 
tion  paralleling  the  Rocky  Mountains  was  considered.  The  location  of 
this  precipitation  band  was  further  east  and  the  amount  of  precipita¬ 
tion  was  greater  than  would  have  been  expected  on  the  basis  of  the 
mean  spring  precipitation  map  alone  (see  Fig.  7).  These  aspects  of 
the  precipitation  band,  therefore,  were  attributed  to  the  effects  of 
the  low  level  trough  extending  northwestward  from  the  North  Dakota 
low.  The  precipitation  was  maintained  by  upslope  flow  and  persisted 
through  April  26,  1969. 

Figure  25  shows  the  precipitation  pattern  generated  after 


six  hours  of  model  time.  The  simulated  trough  had  the  desired  effect. 
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Table  4 


Low-level  Wind  Profiles  -  April  26,  1969 


261200 

GMT 

270000 

GMT 

Stat i on 

Height 

Wi  nd 

He i ght 

Wi  nd 

(gpm) 

(deg/mps) 

(gpm) 

(deg/mps) 

The  Pas 

273 

040/08 

273 

040/08 

591 

050/1 1 

621 

045/08 

101 A 

054/13 

1048 

051/09 

1456 

052/10 

1493 

056/09 

Fort  Smith 

203 

070/04 

203 

120/06 

260 

070/04 

245 

120/06 

659 

118/05 

658 

121/05 

1085 

122/05 

Glasgow 

885 

358/14 

1009 

018/14 

1343 

023/22 

1461 

022/16 

1826 

054/22 

1930 

034/17 

Edmonton 

766 

315/04 

766 

338/02 

1053 

335/11 

1061 

333/03 

1503 

337/10 

1516 

329/04 

1978 

340/10 

1993 

326/05 
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Fig.  2k.  Precipitation  (inches)  pattern,  April  26,  1969. 
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Simulated  precipitation  (inches)  pattern,  April  26,  19&9* 
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The  model  successfully  reproduced  the  band  of  precipitation  leaving 
a  dry  tongue  over  the  foothills,  and  the  amount  of  the  precipitation 
was  enhanced  over  that  produced  in  the  simulation  of  'typical*  late 
spring  upslope  weather.  The  inclusion  of  snow  cover  on  the  mountains 
along  the  western  edge  of  the  grid  and  the  subsequent  lack  of  strong 
upward  heat  fluxes  in  this  area  also  contributed  to  the  production  of 
this  dry  tongue.  However,  the  location  and  quantity  of  precipitation 
was  mainly  due  to  upslope  flow  east  of  the  trough  and  downslope  flow 
west  of  the  trough. 

Figure  26  illustrates  the  simulated  vertical  velocity  pattern. 
The  most  notable  feature  was  the  pattern  of  vertical  velocities  over 

l 

Alberta.  Here,  the  vertical  velocities  were  larger  than  those  found 
in  the  earlier  simulation.  This  effect  was  due  to  strong  low-level 
convergence  along  the  trough  and  to  the  slightly  stronger  low-level 
flow.  Spot  value  comparisons  of  the  resultant  height  of  the  inversion 
base  and  mixed-layer  winds  are  given  in  Table  5.  The  comparisons, 
which  were  severely  limited  by  the  number  of  upper  air  stations  in 
the  area  of  interest  pointed  to  a  successful  reproduction  of  the 
inversion  base  height  and  low-level  flow  patterns  that  actually 
occurred  on  April  26,  1969,  however,  the  data  were  too  limited  to 
draw  any  firm  conclusions. 

Based  on  the  verification  of  the  precipitation  pattern  and 
the  spot  value  comparisons  of  the  height  of  the  inversion  base  and 
low-level  flow  pattern  it  was  concluded  that  the  model  and  the  single¬ 
layer  concept  of  the  PBL  could  adequately  simulate  some  aspects  of  an 
observed  bounda ry- 1 ayer  weather  pattern. 
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Fig.  26.  Simulated  vertical  velocity  (cm  sec"1)  pattern,  April  26,  1969. 
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Table  5 

Inversion  Base  Height  and  Mixed-layer  Winds 


Locat i on 


h* (gpm)  h**(m)  Wi nd*** (deg/mps) 


Wi nds** (deg/mps ) 


Edmonton 

1984 

1965 

332/06 

347/05 

Glasgow 

1881 

1713 

025/18 

018/06 

The  Pas 

1475 

1376 

049/10 

049/09 

^Height  of  the  inversion  base  at  261200  GMT  plus  the 
height  of  the  inversion  base  at  270000  GMT  divided  by  two. 

**Two  grid  point  average  in  the  vicinity  of  the  upper 
ai r  station. 

***Ver tical  mean  of  the  low-level  winds  at  261200  GMT 
plus  the  vertical  mean  of  the  low-level  winds  at  270000  GMT 
divided  by  two. 


CHAPTER  5 


SUMMARY  AND  CONCLUSIONS 

A  three-layer  mesoscale  numerical  model  has  been  used  to 
simulate  boundary  layer  weather  patterns  including  limited  convective 
precipitation  induced  by  a  moderately  moist  east-northeasterly 
circulation  of  arctic  air  over  the  relatively  warm,  sloping  Canadian 
Western  Plains  during  late  spring.  Time-dependent  calculations  were 
limited  to  a  single  'mixed'  layer  by  parameterizing  the  interactions 
of  this  layer  with  both  the  underlying  and  overlying  layers.  Limited 
convective  precipitation  and  the  release  of  latent  heat  within  the 
layer  have  also  been  parameterized.  A  47.6  km  by  47.6  km  grid  mesh 
of  1369  points  covering  the  Canadian  Prairie  Provinces  was  used  to 
represent  the  variables  in  the  mesoscale.  The  governing  equations 
were  solved  numerically  with  surface  forcing  factors  allowed  to 
perturb  the  'mixed*  layer  from  its  initial  conditions  until  resultant 
mesoscale  boundary  layer  weather  patterns  evolved. 

Experiments  using  conditions  'typical'  of  late  spring 
upslope  weather  revealed  that: 

(1)  the  model  atmosphere  provided  a  reasonable  representa¬ 
tion  of  the  atmospheric  structure  associated  with  arctic  outbreaks 
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over  the  Canadian  Western  Plains  during  late  spring. 


(2)  the  mean  spring  topographic  precipitation  pattern  could 
be  reproduced  by  the  simulated  late  spring  upslope  flow  with  limited 
convective  precipitation.  The  'mixed'  layer  representation  of  the 
planetary  boundary-layer  allowed  this  to  be  done  relatively  simply 
and  with  modest  computer  requirements.  Approximately  6.5  minutes  of 
Central  Processing  Unit  (CPU)  time  and  820  page-minutes  of  CPU  storage 
on  an  I BM-36O  computer  were  required  to  simulate  6  hours  of  model  time. 

(3)  the  limited  convective  precipitation  produced  by  an 
east-northeasterly  circulation  over  the  Western  Plains  was  intimately 
related  to  the  configuration  of  the  relief  and  that  upslope  weather 
conditions  persisted  as  long  as  the  low-level  flow  remained  essentiall 
unchanged . 

{k)  late  spring  upslope  weather  is  a  contributing  factor  to 
the  upward  swing  in  precipitation  amounts  over  the  Western  Plains 
during  the  spring  period. 

(5)  mesoscale  PBL  weather  patterns  exert  a  dominant  control 
over  the  location  and  intensity  of  perturbations  in  the  spring  pre- 
ci pi tation  pattern. 

The  sensitivity  of  late  spring  upslope  weather  on  the 
Canadian  Western  Plains  to  surface  heating  has  been  explored.  It 
was  found  that  the  area  and  intensity  of  limited  convective  precipita¬ 
tion  was  significantly  reduced  when  upward  heat  fluxes  were  neglected. 


* 
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that : 


A  simulation  based  on  observed  initial  conditions  showed 


(1)  the  model  was  capable  of  reproducing  a  limited  convective 


precipitation  pattern  maintained  by  upslope  flow. 


(2)  a  low-level  trough  exerts  a  marked  influence  on  the 


location  and  intensity  of  the  precipitation. 


This  numerical  experiment  has  provided  a  preliminary  look  at 
late  spring  upslope  weather  on  the  Canadian  Western  Plains.  However, 
a  number  of  problem  areas  remain.  Parameterization,  initialization, 
and  verification  could  all  be  improved  by  the  incorporation  of 
information  on  a  smaller  scale.  The  lack  of  detail  available  for 
the  mesoscale  distributions  of  drag  and  transfer  coefficients  is  a 
particularly  acute  problem.  Many  more  studies  are  needed  to  complete 
our  understanding  of  late  spring  upslope  weather.  Nevertheless,  it 
is  hoped  that  a  better  understanding  of  this  significant  weather 
phenomenon  has  been  achieved.  More  generally  it  is  hoped  that  this 
simulation  has  provided  further  insight  into  mesoscale  planetary 
boundary-layer  weather  patterns  induced  by  surface  forcing  factors. 
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APPENDIX 


BUILDING  THE  MODEL  ATMOSPHERE 

The  model  atmosphere  was  assumed  to  consist  of  three  layers 
defined  by  distinctive  lapse  rates  of  potential  temperature  (see  Fig. 
3).  The  vertical  structure  of  the  model  atmosphere  was  built  up  from 
surface  parameters  in  the  following  manner: 

(a)  Calculation  of  temperature  T  and  dewpoint  temperature  Td  at 
the  various  levels  in  the  vertical 

(1)  Level  Z 

s 

T  =  T  -  [y  (z  -  Z  )] 

S  O  L  1  5  5  o 

3T 

where  y  =  -  —  was  the  superad  i abat  i  c  lapse  rate  of 
5  o  Z 

temperature. 

Tds  ■  Tdo  -  [5(Zs  •  Zo)] 

3Td 

where  £  =  -  —  was  the  lapse  rate  of  dewpoint  temperature. 
oZ 

(2)  Levels  Z  and  h,  where  Z  =  (Z  +  h)/2 

m  ms 

T  =  Ts  -  [y (Z  -  Zs)] 

where  y  was  the  dry  adiabatic  lapse  rate  of  temperature. 
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Td  =  Tds  -  [5(2  -  ZQ)] 


(b)  Calculation  of  equilibrium  vapour  pressures  eQ,  e^,  e^ 
(dynes  cm”2)  for  levels  Zq,  and  h  using  the  C 1 aus i us-Cl apeyron 
equation 


=  6.11  x  103  EXP 


L.  (J_ 

Rw  273 


’Td* 


where  L  was  the  latent  heat  of  evaporation  assumed  constant  at  the  OC 

value  of  597-3  cal  gm"1.  R  was  the  specific  gas  constant  for  water 

w 

vapour. 


(c)  Calculation  of  the  cloud  base  height  or  LCL 


T  -  Td 

o  o 

y  -  0.  1  7y 


(Pettersen,  1956) 


where  Z  ,  T  and  Td  were  the  initial  values  of  terrain  height,  surface 
o  o  o 

temperature,  and  dewpoint  temperature. 

(d)  Calculation  of  the  vertical  mean  of  the  horizontal  pressure 
grad i ent 

Pressure  values  at  the  various  levels  were  needed  to  convert 
the  temperature  and  vapour  pressure  values  to  potential  temperature 
and  specific  humidity  values.  The  vertical  mean  of  the  horizontal 
pressure  gradient  was  obtained  in  terms  of  the  horizontal  wind  which 
has  been  assumed  constant  with  height.  The  somewhat  arbitrary 
condition  of  balance  between  pressure  gradient,  Coriolis  and  friction 
forces  was  assumed  at  t  =  0.  Vertical  homogeneity  of  the  horizontal 


. 


wind  implies  that  one  pressure  gradient  applies  at  all  levels.  The 
pressure  gradient  is  calculated  by  averaging  the  balanced  equation  of 
motion  through  the  mixed  layer,  while  neglecting  the  flux  of  momentum 
through  the  inversion  surface  which  tops  this  homogeneous  layer. 

Thus,  the  equation  of  motion  becomes  the  Ekman  balance  equation: 


0  =  -  k  x 


fV  -  aVp  +  a 


With  the  assumptions  outlined  above  and  with  ts  =  C^p|v|v  the  above 
equation  becomes: 


aVp 


k  x  f V  - 


Cp|v|v 

(h-zs) 


The  pressure  values  at  the  various  levels  were  then  obtained  in  the 

following  manner.  aVp  calculated  above  was  applied  at  level  Z  . 

m 

Since  density  varies  only  slightly  with  pressure  and  an  explicit 

value  of  pressure  is  required  to  integrate  the  above,  it  seemed  quite 

sufficient  to  use  the  standard  pressure  height  relationship  based  on 

the  U.S.  Standard  Atmosphere  to  obtain  first-guess  pressure  and 

density  fields  at  level  Z  . 

m 


Pq  =  1013.25 


,  6.5  x  10"5Z  *1 

1  -  _ o 

288.16 


6.5  x  10‘5R 


dynes  cm 


-2 


p  =  p 
s  o 


1  - 


Y  (Z  -  Z  ) 

s  s  o 


dynes  cm 


“2 


p  =  p 
m  s 


1  - 


m  s 


_2_ 

YR 


dynes  cm 


-2 


- 
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and  from  the  Equation  of  State 


a  =  a  =  RT  /  p 
m  m  m 


The  pressure  field  at  level  was  then  given  by 


+  IE-  «x  +  IE- 


p  =P  +  -r1-  6x  +  6 y 

rm  rmi  3x  8y  1 


where  p  is  the  (1,1)  grid  point  value  of  the  first-guess  pressure 
field.  Because  the  above  procedure  gave  pressures  somewhat  larger  than 
expected  on  the  basis  of  experience  with  the  type  of  synoptic  situation 
being  considered,  the  value  of  p  was  reduced  by  about  1  percent 
before  the  integration  was  carried  out.  Pressure  fields  at  levels 
Zq,  Z^  and  h  were  then  calculated  using  the  integrated  form  of  the 
hydrostatic  equation. 


where 


where 


p  =  p  EXP 
rs  rm 


RT. 


(Z  -  Z  ) 
m  s 


T.  =  (T  +  T  )/2 
1  s  m 


p  =  p  EXP 
o  rs 


RT 


(Z  -  Z  ) 
s  o 


T2  “  (To+Ts)/2 


p,  =  p  EXP 
rh  rm 


RT 


<zm  * 
m 


T  =  (T  +  T  )/2 
3  m  n 


where 
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(e)  Calculation  of  potential  temperatures  at  levels 

K 


0 


T 


Z  and  h 
s 


where  p^  is  the  reference  level  of  1000  mb  and  k  =  R/C 
to  the  specific  heat  at  constant  pressure  for  dry  air: 


P 


wi  th 


C  equal 
P 


0.2^+0  cal  gm”1  K-1 . 


(f)  Calculation  of  specific  humidity  at  levels  Z  ,  Z  and  h 

o  s 

0. 622e 

q  p  -  0.378e 


(g)  Calculation  of  mean  mixed  layer 
initial i ze  the  model 


to  h)  values  used  to 


6  =  (0  +  0.  )/2 

s  h 

and 

q  =  (qs  +  qh)/2  . 


